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Abstract 

The fluctuations of the jamming coverage upon Random Sequential Adsorption 
(RSA) are studied using both analytical and numerical techniques. Our main re- 
sult shows that these fluctuations (characterized by ctqj) decay with the lattice size 
according to the power-law uqj oc L~^l v . The exponent v depends on the dimen- 
sionality D of the substrate and the fractal dimension of the set where the RSA 
process actually takes place (df) according to v = 2 D-d s • This theoretical result is 
confirmed by means of extensive numerical simulations applied to the RSA of dimers 
on homogeneous and stochastic fractal substrates. Furthermore, our predictions are 
in excellent agreement with different previous numerical results. 

It is also shown that, studying correlated stochastic processes, one can define 
various fluctuating quantities designed to capture either the underlying physics of 
individual processes or that of the whole system. So, subtle differences in the 
definitions may lead to dramatically different physical interpretations of the results. 
Here, this statement is demonstrated for the case of RSA of dimers on binary alloys. 



1 Introduction 



Physical and chemical properties of adsorbed monolayers are being studied with increasing 
interest because their understanding is essential for the rationalization of many phenom- 
ena and processes occurring on surfaces and interfaces, such as adsorption, desorption, 
catalysis, corrosion, wetting, adhesion, diffusion, etc. The equilibrium behavior of such 
overlayers can be described by a Gibbs measure parametrized by the coverage 9 and the 
temperature T. Within this context, the critical behavior of adsorbed films has exten- 
sively been studied [1, 2, 3]. On the other hand, numerous physical processes can be 
modeled as the sequential, irreversible filling of a surface by atoms or molecules. Some 
examples are the reaction at specific sites on a polymer chain, adsorption onto surfaces, 
reaction between groups on adjacent surface sites, etc. [4, 5, 6]. 

Considering the irreversible deposition of particles on a surface, one has two charac- 
teristic time scales: the time between depositions, and the diffusion time of the particles 
on the surface. For very strong interaction between particles and the substrate (chemical 
adsorption), diffusion becomes irrelevant and the venerated Random Sequential Adsorp- 
tion (RSA) model provides an excellent description of the underlying processes (for a 
review on RSA models see e.g. [6]). Under these conditions the system evolves rapidly 
toward far-from equilibrium conditions and the dynamics becomes essentially dominated 
by geometrical exclusion effects between particles. This kind of effects has been observed 
in numerous experiments [7]. 

When RSA involves adsorption on single sites, the case is termed 'monomer filling'. 
Also, processes involving adjacent pairs of sites are referred as 'dimer filling', while ad- 
sorption on larger ensembles of sites corresponds to 'animal filling' [6]. A quantity of 
central interest for the understanding of RSA processes is the asymptotic value of the 
fraction 9 of the total surface occupied by adsorbed objects, which is called the jamming 
coverage 9j = 9{t — > oo). Within this context, the RSA of needles (or linear segments) on 
homogeneous, two-dimensional samples, has very recently attracted considerable interest 
[8, 9]. Particular attention has been drawn to the interplay between jamming and per- 
colation [8, 9, 10]. Of course, the percolation problem is a topic of enormous interest by 
itself, due to their applications not only in statistical physics but also in many other areas 
such as the study of disordered media, fluids in porous materials, systems of biological 
and ecological interest, etc. [11, 12, 13]. , a great progress in the field of the statistical 
physics of far-from equilibrium processes could be achieved by establishing links between 
RSA and percolation [8, 9, 10]. 

Percolation is essentially a geometrical critical phenomena. The percolation transition 
is related to the probability of occurrence of an infinite connectivity between randomly 
deposited objects, as a function of the fraction p of the substrate occupied by the objects 
[14]. At the critical point (p c ), the percolation cluster in d— dimensions, is a random fractal 
of dimension Dp = d — fl/v, where f3 and v are the order parameter and correlation length 
critical exponents, respectively. Close to criticality, the probability P to find a percolating 
cluster, on a finite sample of side L can be described by means of an error function [13] 




(1) 



where A is the width of the transition region. It is well known that the width vanishes in 
the thermodynamic limit according to [13] 



AocL~. (2) 

Eq. (2) is very useful because it allows for the measure of the correlation length exponent 
v that governs the divergence of the correlation length as £ oc \p — p c \~ u (for random 
percolation with d = 2 one has v = 4/3). 

Very recently it has been suggested that the jamming probability and the fluctua- 
tions of the jamming coverage may obey relationships similar to Eqs. (1) and (2) [8], 
respectively. Also, for the jamming upon RSA of needles in two dimensions the value 
vj = 1.0 ± 0.1 has been reported [8] and this figure is independent of the aspect ratio 
of the needles. Furthermore, early numerical results of Nakamura for the RSA of square 
blocks are also consistent with vj ~ 1 [15], while Kondrat and Pekalski [9] have reported 
vj = 1.00 ± 0.05 for the RSA of segments on the square lattice. Since the obtained values 
for the exponent are independent (within error bars) of: i) the length of the segments (for 
all a = 1,2,. ...,45) [9], ii) the aspect ratio of the needles [8] and iii) the size of the square 
blocks [15], it has been suggested that vj is a good candidate for an universal quantity of 
the jamming process [9]. 

The aim of this manuscript is to provide a qualitative derivation of Eq. (2) for the 
case of RSA on heterogeneous media. It is shown that the exponent vj can be obtained 
as a function of the dimensionality D of the space and the fractal dimension df of the 
subset of sites where the RSA process actually takes place. Our main result 

vj = 2ihdi (3) 

provides a solid ground to previous numerical data [8, 9, 15]. In fact, these data were 
obtained in D = 2 and df = 2, so it follows straightforwardly from Eq. (3) that vj — 1 
exactly. Furthermore, in this work, the validity of the proposed relationship (3) is verified 
by means of extensive numerical simulations, using both homogeneous substrates as well 
as different random fractals. 



2 Definitions and Approaches 

When studying RSA on homogeneous samples the definition of the jamming coverage 
and its fluctuations is straightforward, since one has to deal with a single stochastic 
process. However, RSA on nonhomogeneous substrates may involve the treatment of at 
least two correlated stochastic processes: the selection of the particular substrate where 
the deposition is going to take place, and the RSA process itself. In relation with the 
first process we will define the substrate system as the set {^a}, composed of M different 
independent substrates labeled by the index A = 1, M. We consider that n independent 
RSA processes are performed for each element A\ of the substrate system until the jammed 
state is reached. We call B^ (k = 1, ...,n) the n different configurations adopted by the 
entities adsorbed on top of the substrate. The set {B^} taken for the values of k (for 
each and all substrate) will be referred to as the RSA system. 



The jamming coverage 9j is a relevant intensive quantity that takes the value 9^ 
when evaluated over the configuration . Let us explicitly consider two contributions 
to 9 j given by 



9 J — d{sub} + $9{Rsa} , 



(4) 



It is also assumed that the first term gives the most important contribution to 9j, that 
can be evaluated according to 
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Then, the fluctuations of this term are given by 
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where is the average value of 9j, taken over all measurements and configurations 
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The second term in Eq. (4) has a zero average for a fixed substrate. Some fluctuations 
in 9 must also appear from this term and can be characterized by the following average 
[16] 
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af^ being given by 
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It should be stressed that <7{ sub } and <7{R sa } account for the fluctuations on the substrate 
and RSA system, respectively. So, it is expected that these quantities will describe the 
relevant physical behavior related with both sources of randomness. 

Finally, the total root mean square deviation (RMS) a$, that is expected to describe 
the fluctuation of the whole system, is given by 
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In general, different RMS's are related. In the simplest case when both contributions 
in Eq. (4) are statistically independent, the total fluctuation of 9 is given by 



a ° = \J a {Sub} +(T {Rsa}- (H) 

As will be shown later, Eq. (11) is essential for the analysis of numerical data obtained 
using finite systems and their corresponding extrapolation to the thermodynamic limit, 
because each term may behave according to different laws. Therefore, simply measur- 
ing erg, one may obtain inaccurate results, undesired crossover effects and consequently 
misleading physical interpretations. 

Let us now evaluate the fluctuations associated with the RSA system using the defini- 
tion given by Eq. (8). Let = 0, 1 define the occupational state of the site i. N = Yli=i n i 
is the number of particles deposited in a volume Q = L D of a D-dimensional lattice with 
periodic boundary conditions. The two-point correlation function for a homogeneous 
substrate is given by 

G(i,j) = (n i n j )-(n i )(n j ) (12) 

where i,j labels sites and () means average taken over the RSA system in which N varies. 
There is a fundamental relationship between G and the fluctuation in the number of 
particles given by 

L D 

o%=Y,G(i,j). (13) 

*j=i 

For the case of a unique homogeneous substrate it is well known that, when the 
correlation length is much shorter than the size of the system, one has 

a 2 N = gN A . (14) 

Here g = ~Y^j=i G(i,j) is independent of Na, the number of sites were RSA takes place. 
The fluctuation of the density (9) in the thermodynamic limit can be obtained after 
dividing both sides of Eq. (14) by a L 2D so that 

a {Rsa} oc L~%, (15) 

Let us now generalize this result in order to show so that it still holds for substrates 
having sites where deposition is forbidden. The existence of these blocked sites on the 
adsorbing surface breaks the translation symmetry of the substrate, which was a necessary 
condition in the above deduction. However, an equation analogous to Eq. (15) can still 
be found if we restrict ourselves to those cases where although each substrate A\ is 
nonhomogeneous, the set {A\} is homogeneous as a whole, i.e., there is no preferred sites 
on the average. Stated in a more quantitative way, the following calculations are valid 
for systems where the average of n,i over substrates is independent of %. Then, using a 
procedure analogous to that employed in the derivation of Eq. (8), one can generalize Eq. 
(13) by writing 

A=l i,j=l i,j=l A=l 



where the expression within parenthesis could be a measure of correlations over the sub- 
strate system. However, for a fixed site % there may exist substrates in the adsorptive 
matrix where the deposition is forbidden, so that G\{i,j) = for all j. Then, since not 
all the M substrates would be contributing in the sum over A, simply dividing by M one 
cannot obtain the real average of G over A. In order to have a proper average of the RSA 
correlation function, it is necessary to introduce a factor X i} such that XjM is the actual 
number of substrates with G\(i,j) ^ 0. It should be noticed that for systems that are 
homogeneous as a whole one has Xi = X, so that 

i,j=l A=l ij'=l 

In order to evaluate X, the fraction of substrates where a site is not blocked, one has 
to count the number of non-zero terms in Eq. (16). Adding first over j one obtains 

L D 

Notice that this procedure implies that g\(i) = for blocked i-sites in a fixed substrate 
A, while g\(i) ^ otherwise. The remaining double sum 

M,L D 

E 9x(i), (19) 

A=l,i=l 

can be computed fixing A and running % over the sites. Then one obtains that the number 
of contributing terms is exactly the number of active sites in this substrate. Therefore, 
the summation over A has ML d f terms, where L d f is the average number of active sites 
of the lattice where RSA actually takes place. 

When % is fixed and the sum given by Eq. (19) runs over A, it yields MXi terms 
(instead of M). Then the summation ^2f =1 MXi has ML D Xi terms. Comparing now 
both results and using again the fact that there is no preferred sites in the global system, 
it follows that X = -jp. It should be noticed that this correction factor is not totally 
unexpected. In fact, for a substrate system that is homogeneous as a whole, one expects 
that the fraction of substrates with a fixed site i active Xj, should equal the average 
proportion of active sites given by L d f / L D . 

Now, replacing X in Eq. (17), it follows 

If the substrate system is homogeneous as a whole, then G* should be invariant under 
translations, because the average of G over A should have the same symmetries than 
those of the substrate system. Thus, g*(i) = Y^=i G*(hj) m ust be site-independent (i.e., 
g*{i) = <7o)j an d then Eq. (20) becomes 




(21) 



Finally, if the correlation length ^ Rsa ^ associated with G* for the RSA process is short 
enough (^ Rsa y << L), g$ is L-independent. Then the fluctuation of the density (9) can 
be obtained from Eq. (21) dividing by L 2D , so that 



It should be stressed that Eqs. (22) and (23) are quite general relationships valid for 
substrate systems that are homogeneous as a whole. Furthermore, the condition that the 
correlation length of the RSA process should be smaller than the system size is usually 
valid for jammed states, where the correlation length is very short. 

Eqs. (22) and (23) also holds for homogeneous substrates, with L d f = L D and G\ = 
G, so v = j^. These equations are also valid for nonhomogeneous random substrates 
where the dimensionality D of the space may be different to the dimensionality of the 
subset of sites where the RSA process actually takes place. It is also very interesting to 
notice that, using these relationships it may be possible to evaluate df performing RSA 
both in numerical simulations and actual experiments. Furthermore, existing numerical 
simulations performed in D = 2 dimensions with df = 1 are in excellent agreement with 
Eqs. (22) and (23) (notice that for these conditions it follows straightforwardly from Eq. 
(23) that v — 1 exactly): i) Nakamura [15] has reported vj ~ 1 (RSA of square blocks), 
ii) Vandewalle [8] have reported vj = 1.0 ± 0.1 (RSA of needles), and iii) Kondrat and 
Pekalski [9] have reported vj = 1.00 ± 0.05 (RSA of segments). 

In addition to these promising results, we will provide more astringent tests of the va- 
lidity of Eq. (23) performing numerical simulations of the RSA of dimers on homogeneous 
substrates in D — 1, 2 dimensions and as well as using fractal substrates. 

3 Details on the Numerical Simulation of RSA of 
Dimers 

3.1 RSA of dimers on Binary alloys 

The first set of simulations are performed for the RSA of dimers (the RSA system being 
the adsorbed atoms) on a binary alloy (BA) annealed at a given temperature (substrate 
system). The BA is simulated using the isomorphism with the Ising model [17], namely 
spin-up = A-species and spin-down = Z?-species. The square lattice of side L in D — 1, 2 
dimensions with nearest-neighbor (NN) interactions will be considered. The BA is in 
contact with a thermal bath at temperature T. The system is assumed to obey Kawasaki 
dynamics [17], so that the density of species A and B is conserved with pa = Pb = 1/2. 
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The Hamiltonian (H) is given by 

H = E - Jj2 s i s i ( 24 ) 

(id) 

where E = N(e AA + e BB + 26 A b), J = — (l/4)(e^ + ess - 2e AB ), e X y is the interaction 
energy between X- and F-species, and Sj = ±1 indicates A, B sites so that ^ s i — 
[18]. The Monte Carlo time step (MCS) involves L D trials, such that each species of the 
sample is selected once on average. Stationary configurations of the BA are obtained 
after disregarding 10 5 MCS. It is well known that for D — 1 the system is not critical, 
exhibiting a (trivially) ordered phase only for T = 0. However, in D = 2 dimensions the 
system undergoes an order-disorder transition at T c = 2.269..., where T c is the Onsager 
critical temperature (the temperature T is measured in units of J setting Boltzmann 
constant at unity). Below this transition temperature, the dynamics of the underlying 
Ising model can be very slow because the system gets trapped on different metastable 
states separated between each other by high energy barriers. In order to have meaningful 
and systematic results, we have taken the same initial configuration for the simulation of 
each low temperature substrate (T < T c ). The starting configuration chosen was the one 
with lowest energy, in which A- and 5-species are segregated into two identical domains, 
separated by non-defective (straight) domain-walls. 

RSA of dimers on the BA is assumed to take place on top of NN sites occupied by 
unlike species, i.e., AB-sites, adsorption on A A- and BB-sites being forbidden. Also 
multiple occupation of AB-sites is forbidden. These assumptions are based on the disso- 
ciative chemisorption of diatomic molecules (O2, H 2 , N 2 , etc.) on binary catalysts. It is 
well known that in the absence of dimer's diffusion there is a jammed state [19], so that 
the relevant quantity is the jamming coverage (0j). In the case of homogeneous (D = 1) 
lattice one has the exact result 9j = 1 — e~ 2 ~ 0.86466472 [6]. Also, in the homogeneous 
(D = 2) lattice one has 0j ~ 0.906 [6]. For the present case of a BA one expects the 
coverage to depend on the temperature at which the substrate has been annealed. Also, 
it should satisfy 9^ B ^ < 6j, due to the additional constraint that dimmers can only be 
adsorbed on specific pairs of sites of the lattice. 

3.2 RSA of dimers on fractal surfaces 

We have also studied the RSA of dimers on stochastic fractals such as the diffusion front 
[11, 12] and Ising clusters [17]. It is important to stress that these chosen fractals satisfy 
the required property of invariance due to their stochastic nature, so that Eqs. (22) and 
(23) are expected to hold. 

Considering a BA at T = T C) the greatest cluster made up taking NN sites occupied 
by the same species, is selected for different lattice sizes. This kind of substrate, called 
spin clusters of the Ising model (SCIM), are fractals with a fractal dimension given by 

df = D - (25) 

where D is the Euclidean dimension, (3 and v are the order-parameter and the correlation- 
length critical exponents, respectively. For D = 2 one has f3 — | and v — 1, so that 
d f = 15/8. 



In order to simulate the RSA of dimers on SCIM's, a site of the cluster is selected at 
random. If that site is empty, a NN site is also selected at random. The adsorption trial 
is successful only if that second selected site is also empty, otherwise if that site is either 
occupied or lies outside the fractal, the trial is disregarded. Therefore, double occupation 
of sites belonging to the fractal is forbidden. 

RSA along the perimeter of SCIM's has also been studied. For this purpose, only 
adsorption events of dimers taking place on two NN sites, such us one of then belongs to 
the cluster and the remaining one is outside it, are considered. 

Another type of stochastic fractal used for RSA simulations is that generated by a 
diffusion front. In fact, it is well known that the properties of the diffusion front [21, 22, 
23, 24, 25] are closely related to those of the incipient percolation cluster [11, 12, 13]. 
The front is a (stochastic) self-similar fractal of dimension [20] df = 10/7 ~ 1.42857 
[21, 22, 23]. 

In order to obtain diffusion fronts suitable for RSA of dimers, we have simulated the 
diffusion of particles at random. Hard-core interactions, on a D = 2 square lattice of size 
L x L, has been considered. There is a source of particles at the first row of the lattice 
y — 1,1 < x < L kept at concentration p(y) = 1. Also, at the last row y = L, 1 < x < L, 
there is a well such as p(L) = 0. So, there is a concentration gradient along the source- 
well direction, while along the perpendicular ^-direction periodic boundary conditions are 
imposed. As the concentration p(y) of particles depends on the position, decreasing from 
the source to the well, one actually has a gradient percolation system. The structure 
of the diffusion front is identical to the structure of the hull of the incipient percolation 
cluster [21, 26]. Furthermore, the concentration of particles at the mean front position yf 
is the same as the percolation threshold p c , so that p{yf) = p c [21, 27]. 

Considering that particles are connected to their first neighbors while empty sites 
are connected to both their first and second neighbors, the system can conveniently be 
described by means of the following geographical analogy: The land is the set of particles 
connected to the source and the sea is the set of connected empty sites not surrounded by 
land. Also, an island is formed by groups of particles not connected with the land, while 
lakes are formed by connected empty sites surrounded by land. Finally, the seashore 
is the part of the land in contact with the sea that can be identified with the diffusion 
front, which is a self-similar fractal. 

RSA of dimers on diffusion fronts has been simulated using two rules: Rule I is that 
used for the case of RSA on the SCIM perimeter (see above); while using Rule II one 
only allows the adsorption on sites of the diffusion front, disregarding adsorption trials 
on already occupied sites of the front and sites outside the fractal. 

4 Simulation results and discussion 
4.1 RSA of dimers on Binary Alloys 

The simplest test for Eq. (23) corresponds to the case D — df — 1. Since for BA we 
are only allowing deposition on unlike sites, this example can be realized taking a one di- 
mensional BA with nearest-neighbor repulsion between unlike species (antiferromagnetic 
Ising system). In fact, considering the ground state at T = 0, corresponding to a fully 




Figure 1: Log-log plots of <J{R sa } versus L obtained the case of RSA on BA's in 1 dimension 
(D = 1). BA's at two different temperatures, T = and T = oo, are considered as shown 
in the figure. 

ordered sample, one has that the BA is irrelevant and the system is equivalent to the 
standard RSA of dimers in D — 1. On the other hand, for T > one has disordered 
samples with an homogeneous distribution of blocked sites, i.e., AA- and 55-sites where 
dimers are not adsorbed. Particularly, at T = oo one has a fully disordered substrate. 
The obtained results are shown in Fig. 1 as log-log plots of <7{R sa } versus L. The best fits 
of the data give v = 2.00 ± 0.01 and v = 2.01 ± 0.01 for the cases T = and T = oo, 
respectively. Both figures are in excellent agreement with the prediction of Eq. (23) for 
D = df = 1, namely v — 2. 

The second example considered corresponds to RSA on a D = 2 BA with NN attrac- 
tive interactions between alike species, i.e the ferromagnetic version of the Ising model 
with conserved order parameter. In this case one can also tests the validity of Eq. (11) by 
measuring all the involved terms according to Eqs. (6), (8) and (10), respectively. Fig. 2 

shows log-log plots of a,? versus \J<?\ Sub } + °"{i? sa } obtained for the case of RSA of dimers 

on two dimensional BA's. Data was taken at different temperatures and using lattices 
of different size. The excellent straight line obtained, corresponding to 150 independent 
measurements, strongly supports the validity of Eq. (11). 

Also, Fig. 3 shows plots of o"e{_Rsa} versus L for the RSA of dimers on BA's obtained at 
different temperatures. It is found that the power law decay predicted by Eq. (22) always 
holds allowing us to determine the exponents v. At low temperatures, below T c , one has 
that the BA segregates into domains of alike species containing a certain density of unlike 
species (impurities) trapped into the bulk, which increases when the temperature is raised. 
Since dimer adsorption on the bulk of the domains is not possible due to the adsorption 
rule that have been imposed, at very low temperatures the RSA process is essentially 
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Figure 2: Log-log plots of a e versus y o~\ Sub } + °"{_r S(1 }> as suggested by Eq. (11). The dif- 
ferent terms involved in Eq. (11) were obtained according to Eqs. (6) and (8) (horizontal 
axis) and (10) (vertical axis), respectively . The straight line has slope unity. 

restricted to the interface between domains. In this case one has df = 1 and Eq. (23) 
predicts v = 2/3, in excellent agreement with the results obtained fitting the curves, as 
shown in the inset of Fig. 3. Increasing the temperature and particularly close to T c , the 
density of impurities located into the bulk of the domains increases. Consequently, an 
increasing number of pairs of sites becomes available for the RSA of dimers. However, 
not all these sites contributes to the fluctuations of the jamming coverage. The simplest 
example is the case of a single impurity surrounded by unlike species where only one dimer, 
with four possible orientations, can be adsorbed. Therefore, fluctuations of the jamming 
coverage are only relevant when adsorption takes place in rather complex arrangements 
of impurities that are present in the domains close to T c . For these reasons one observes a 
smooth increase of v when approaching the critical point from below, as shown in the inset 
of Fig. 3. Finally, for T > Tq one has that adsorption sites, given by nearest-neighbor 
pairs of unlike species, are homogeneously distributed on the sample with D = 2 and 
df = 2. For this case, Eq. (23) predicts v — 1 in excellent agreement with the numerical 
results shown in the inset of Fig. 3. It should be noticed that the smooth variation of 
v, observed in the inset of Fig. 3 when approaching T c , may be due to finite-size effects 
that hinder the evaluation of the actual exponents. If this is the case, the exponents 
2/3 < is e ff < 1, may be considered as effective size-dependent exponents. 

In order to further support the above discussed interpretation of the evaluated values of 
the exponent z/, numerical simulations using samples having controlled interface roughness 
and concentration of impurities, have been performed. The starting substrate is a ground 
state of the BA (T = 0) in D = 2. The half-left (-right) side of the sample is filled with 
A- (B-) species, so that a perfectly flat interface of A£?-species runs along the middle of 
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Figure 3: Log-log plots of <7{R sa } versus L obtained at different temperatures performing 
in each case n = 500 and M = 500 measurements. The inset shows the temperature 
dependence of the exponent v (temperature is measured in units of J). 

the sample. Since RSA on this substrate has zero rms, the original interface is modified 
introducing kink-like defects with probability p k . The resulting interface where the RSA 
process actually occurs is homogeneous with fractal dimension df — 1, so that according 
to Eq. (23) one should expect an exponent v — 2/3. This prediction is again in excellent 
agreement with the numerical results obtained taking p k = 0.005, 0.05, 0.3 and 0.5, as 
shown in Fig. 4(a). Of course, for p k = 0.005 and small lattices (L < 100), finite size 
effects are observed due to the finite probability of having a perfectly flat interface without 
fluctuations. 

In addition to the test described previously, samples with a fixed interface roughness 
(Pk — 0.3) in the example shown in Fig. 4(b) were decorated with a controlled density 
(pi) of impurities uniformly distributed in the bulk of the domains. Fig. 4(b) shows that 
the presence of impurities not only causes the fluctuations to increase but also v — > 1 for 
large values of pj, as expected for D = 2 and df = 2 according to Eq. (23). The value 
v = 0.70 ± 0.03 obtained for pi — 1/8 is very close to the figure v = 2/3 corresponding to 
pi = since a majority of isolated and small clusters of impurities can not significantly 
influence neither the fluctuations nor the exponent. In fact, only more complex clusters of 
impurities, as those formed for pj = 1/3 and pj = 2/3 in the example shown in Fig. 4(b), 
allow for a large variety of adsorption configurations with an appreciable enhancement of 
the fluctuations of the RSA process that causes a noticeable effect on the exponent. 

4.2 RSA of dimers on fractal surfaces 

Fig. 5 shows log-log plots of cr^ Rsa y versus L obtained upon RSA of dimers on SCIM's. 
In this example one has D = 2 and df = 15/8, so that the value v = 16/17 ~ 0.941176 
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Figure 4: Log-log plots of <T{R sa } versus L obtained for a simple one- dimensional interfacial 
model embedded in a two dimensional lattice. The parameter allow us to choose 
different interface roughness, while pi accounts for the density of two-dimensional defects 
in the bulk. For additional details on the model see the text, (a) Data corresponding to 
pi — and different values of the roughness, as listed in the figure, (b) Data obtained 
keeping the roughness pk = 0.3 fixed and using different values of pi, as listed in the 
figure. 



o DF Rule I, u=0.77 ±0.01 

v DF Rule II, u=0.77 ±0.01 

□ SCIM Bulk, v= 0.941± 0.007 
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Figure 5: Log-log plots of o\r Rsa } versus L for the case of RSA on random fractals. Con- 
sidering Spin Clusters of the Ising model (SCIM) at T c , it is found that Eq. (23) is valid 
when the RSA process is done on both, the bulk and the perimeter cluster. Also, results 
obtained upon RSA on the fractal generated by diffusion fronts (DF) with two different 
adsorption rules, are shown. For details on the SCIM and the adsorption rules see the 
text. 

is expected according to Eq. (23). The results obtained fitting the data corresponding 
to RSA on the bulk of the cluster and on its perimeter are v = 0.941 ± 0.007 and 
v = 0.93 ± 0.02, respectively. In the latter case for small lattices (L < 60), finite size 
effects are observed as in the analogous percolation problem [13]. Also, Fig. 5 shows 
the results for the RSA of dimers on a diffusion front in D = 2. This is an interesting 
example to test the validity of Eq. (23) using an stochastic fractal with a well known 
fractal dimension df = 10/7. For this example one expects v = 0.777... while the best 
fit of the numerical data shown in Fig. 5 gives v = 0.77 ± 0.01, using two different 
adsorption rules as discussed above. These results are also in excellent agreement with 
Eq. (23) giving further numerical support to the analytical relationship derived in Section 
2. 

4.3 Interplay between fluctuations of correlated stochastic pro- 
cesses 

It should be stressed that in order to capture the physics of the RSA process, all the 
fluctuations measured in the previous sections were evaluated with the aid of Eqs. (8) 
and (9). The aim of this subsection is to show that a quite different physical picture can be 
obtained measuring the fluctuations linked with the substrate system. Let us recall that 
for this purpose one also has to perform RSA measurements, however fluctuations have 
to be evaluated using Eqs. (6), which within the context of the present work, captures 



the physics of the underlying substrate where the RSA process takes place. 

As an example, the discussion will be restricted to the case of RSA of dimers on a 
BA in D = 2, in order to perform comparisons with the data obtained in the simulations 
shown in Figs. 2 and 3. Considering a single configuration of the alloy and according to 
the adsorption rules, one has that 9j is essentially a measure of the density of AB-sites 
present in the configuration. Furthermore, such pairs (which are known as the 'broken 
bonds' in Ising spin language), contribute to the internal energy per site of the BA (u). 
Applying the fluctuation-dissipation theorem [28] it is possible to show that the specific 
heat of the BA (Cy = (fjOy) * s gi ven by the fluctuations of u, namely 

L 2 " ( Ux - (u)Y IS 
C v = f22^ M _ x =T2 a u (26) 

A=l 

where averages are taken over the set {^a} of different configurations of the BA. On the 
other hand, for the D = 2 Ising model one has that 

C V (T) oc \T-T c \ a , (27) 

with a = 0, i.e. a logarithmic divergence at criticality. Analogously, one can define an 
RSA 'susceptibility' x f° r the jammed state as 

x ~ L 1^ M-l ~ {Sub} ' { } 

A=l 

where it is expected that, if the system conformed by the deposited particles would follow 
(or 'copy'), in some way, the structure of the underlying substrate, x would follow the 
same behavior that Cy. 

Fig. (6) shows that plots of x versus T, obtained using lattices of different size, exhibit 
clear peaks close to the critical temperature of the underlying B A, resembling the behavior 
of the specific heat in finite samples. In fact, the peaks are shifted and rounded due to 
operation of finite-size effects. This behavior is the typical one for a second-order phase 
transition that implies the existence of a diverging correlation length £ when approaching 
criticality according to [17] 

Z(T)<x\T-T c \- v *, (29) 

where v* — 1 is the correlation length exponent of the Ising model. Using finite-size 
scaling arguments one can set V-I V *\T -T c |«l [17], then replacing into Eq. (27) with 
a — 0, using Eq. (28) and assuming that close to criticality one has x ~ C v , it follows 

Xmax(L) oc ln(L), (30) 

where Xmax is the maximum value of x that can be obtained from the peaks shown in 
Fig. 6. The results shown in Figs. 7(a) and 7(b) confirm the divergences of x expected 
according to Eqs. (27) and (30), respectively. So, this finding shows that there is an 
additional (divergent) correlation length of the RSA process that is associated to the 
substrate, and it can be captured by measuring x- 
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Figure 6: Plots of \ versus T obtained for the RSA of dimers on BA using lattices of 
different size L (T in units of J). The data exhibits the typical behavior characteristic of 
a second-order phase transition. 

From Figs. 3, 6 and 7, as well as according to the above discussion, it follows that 
the RSA of dimers on BA's provides an interesting example of the interplay between two 
correlated processes whose respective fluctuations obey different functions. It is worth 
mentioning that, in spite of exhibiting quite different behavior, both functions capture 
the criticality of the substrate. In fact, the fluctuations of the jammed state, as measured 
according to Eqs. (6) and (28), reflect the critical behavior of the BA through the 
relationship between x an d C v while, on the other hand, using Eqs. (8) and (9) the 
criticality of the BA becomes evident through the jump of the exponent v observed close 
to T c (see Fig. 3). It is also interesting to remark that the relationship that relates the 
different fluctuations of both processes, given by Eq. (11), does not show such a critical 
behavior at all, as shown in Fig. 2. 

In order to further analyze these results, let us define 0* = cr{R sa }l l(J {Sub}, such as 0* 
gives a measure of the relative intensity of the RSA fluctuations as compared to those of 
the substrate. Then, using Eqs. (28), (22) and (30), it follows that 

0* = 0(T)/(L) = 0(T) L (31) 

where the L-dependence behavior of 0* appears explicitly through the function f(L) : 
while 0(T) accounts for the temperature dependence. Fig. 8 shows that log-linear plots 
of 0(T) = 0*//(L) versus T exhibit an acceptable collapsing. At both sides of Tc the 
function 0(T) approximately follows a logarithmic behavior, resembling the temperature 
dependence of \ near T c . Moreover a clear jump appears at T c , due to the change in the 
fractal dimension of the adsorbing set of sites of the substrate (df — 1 for T < T c and 
df = 2 for T > T c ) observed in the finite samples used in the simulations. This result 
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Figure 7: Data corresponding to the RSA of dimers on BA's. (a) Linear-log plot of Xmax 
versus L. (b) Linear-log plot of x versus \T — T c \ (in units of J) obtained using lattices 
of size L = 256. Both plots show that x reflects the divergences of the specified heat of 
the underlying BA. 
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Figure 8: Plot of versus T (T in units of J), obtained using lattices of different size L. 
The jump observed in occurs at T c (dashed line). The solid lines have been drawn to 
guide the eyes. 

implies that the temperature dependence of the structure of the substrate prevails over 
that of the RSA process, and that the fluctuations due to substrate also prevails over 
those due to the RSA process (notice that f(L) — > when L — > oo). 

5 Conclusions 

In this paper the behavior of the fluctuations of the jamming coverage upon RSA process 
on homogeneous and nonhomogeneous substrates has been studied. Pointing the atten- 
tion to the RSA process and applying the definition given by Eq. (8) in order to measure 
the fluctuations, we have shown both analytically and by means of numerical simulations, 
that the jamming fluctuations behaves as ag r oc L^ 1 ^ 1 , where vj is given by vj = 2D 2 _ df , 
D and df being the dimension of the lattice and that of the active sites where adsorption 
actually takes place, respectively. These results are suitable to describe systems charac- 
terized by a short-range correlation length of the RSA process, that may take place on 
both, homogeneous and nonhomogeneous substrate. From our derivation of Eq. (23) it 
follows that vj = 1 for D = 2 and df = 2, in accordance with results published previously 
by various authors [8, 9, 15]. 

It should also be noticed that Eq. (23) allow us to measure the fractal dimension 
of the different adsorption sets where deposition takes place according to the specified 
adsorption rules. Summing up, our results not only point out that a careful treatment 
of the fluctuations of correlated processes is necessary in order to capture the desired 
physical behavior, but also provide a tool for the evaluation of fractal dimensions using 
RSA experiments. 
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